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Q\ ' Based on an unconventional Fermi liquid model, we present several results on the optimally 

doped and overdoped cuprate superconductors. For the normal state, we provide an analytic demon- 
stration, backed by self-consistent Baym-Kadanoff (BK) numerical calculations, of the linear in T 
resistivity and linear in 1/e optical conductivity, provided the interacting Fermi liquid has strong 
peaks in its density of states (van-Hove singularities in 2 dimensions) near the chemical potential 
fx. Moreover, we find that the interactions tend to pin these strong density of states peaks close to 
\x. We show that the low energy dependence of xmmp has a fermionic origin. We obtain particu- 
larly high transition temperatures T c from our BK-Eliashberg scheme by introducing an ansatz for 
the fermionic susceptibility of the carriers. We postulate that the latter is enhanced in an additive 
manner due to the weak antiferromagnetic order of the CuCh planes. We have obtained a d x 2_ y 2 
O ' gap with T c > 120°K for n.n. hopping t = 250meV. 

o 

The nature of the many-body state of the cuprate superconductors is a core question for the understanding of these 
materials |l] [|]. In this Letter, we show that a minimum unconventional Fermi liquid model accounts in a natural, 
comprehensive and internally consistent manner for several normal state characteristics. The introduction of an ansatz 
for the susceptibility of the carriers further allows us to obtain particularly high transition temperatures T c . 
d . Our starting point is the 2-dimensional Hamiltonian 
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H — ^2 £fe C k,a C k^ + — J! C \+q,cr C \>~q,a> C k' ,<?' c k,a + "y C iLg,<r4' -q-<T C K ,-<T C k,<T ■ (1) 

£H , k,a k,k f ,q,(T,a' k,k' ,q,a 

c\ is an electron creation operator and ej, = — 2t{cosk x + cosk y ) — At' cosk x cosk y — 2t" (cos2k x + cos2/cj,). A a is 
the strength of a contact Coulomb interaction, and U the strength of a paramagnon interaction. A small negative 
^vq ' A would mimic an isotropic phononic interaction. 

We have developed a self-consistent many-body treatment of the above Hamiltonian. We write a set of diagrams 
7-H ' within the frame of a Baym-Kadanoff conserving approximation Q] , summing bubble and ladder diagrams Q . There 
is a free energy functional $[G] of the Green's function G, such that the self-energy E is given by the relation 
,— ' E = 5&{G]/5G. We thus obtain a set of self-consistent equations for G(k,e n ) and E(fc, e„) : 

m ■ 

O ' 1 T v-^ 

G(k,e n ) = - ■ — , E(fc,e n ) = --^ V V{q,w m ) G(k-q,e n -u m ) . (2) 

rr\ le n + M - e k - E(fe,e n ) N z *-f 

fi is the chemical potential and the Matsubara frequencies are c n — (2ti -\- l)7rT and uj m — 2uittT for fermions and 
bosons, respectively. This system of equations is solved numerically, similarly to e.g. We work with a given 

number M of Matsubara frequencies and a NxN discretization of the Brillouin zone (M = 256—480 and N > 64). The 

potential V(q,(J m ) = V ex (q,uj m ) - V H (q,u m ) ,V ex (q,u m ) = A {I { { A"+u°y^i{q q ^Z) i-A'ixl(q^ m ) ~ A ° ,V H (q,u m ) = 

i-(A U +u ^xtq^) + l-A^x^q^ ) & Vex(Q> 0Jm ) mciu des vertex corrections due to the internal dressing of bubbles 
with An. Our potential is an extension of the fluctuation-exchange (FLEX) approximation of Bickers, Scalapino and 
White ||, as can easily be seen by taking A a — > 0. Moreover, Xo{q, w m ) = —(T/N 2 ) J2k e + 9> £n + UJ m)G(k, e„) . 

All the convolution operations are done by using the Fast Fourier Transform (FFT), in order to cut down calculation 
time. We use Pade approximants p3| to analytically continue our results to the real frequency axis. 

The numerical solution of the many-body system yields always an interesting result for the self energy at finite 
temperature. 7mE(fc, e) turns out to be essentially linear in energy in an interval t\ < e < 62. However it does have 
the correct parabolic Fermi liquid bevahior for e — * 0. Furthermore, the energy interval of linear behavior expands as 
the energy e v n at which are located the van-Hove singularities (vHs) at the points q = (±7r, 0), (0, ±7r) approaches 
the chemical potential /1. 

We give an analytic derivation of this result. At finite temperature ZmE(fc, e) is given by fl4|| : 

ImT, R (k,e) =J2lmG R (q,e - uj)ImV R (k - q,uj) {coth(w/2T) + tanh((e - w)/2T)} . (3) 
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Taking ImG R (k, e) = — wS(sk, e ), Sk,e = e + M — e k > we obtain 



ImY< R (k,e) = -TrJ2lmV R (k - q,s qte ) {coth{s q , e /2T) + tanh((e 9 - fi)/2T)} . (4) 

Setting ImG(k, e) equal to a delta function is a reasonable approximation for this purpose, since numerically ImT,(k, e) 
is very small compared to the band energy for small couplings. Then, the difference of imX(/c,e), as seen in our 
numerical calculation, for small and large coupling constants is quantitative rather than qualitative. 

We take ImV R (q, x) — X^^Lo Vq 2n+1 \o) x 2n+1 /(2n + 1)!. This is true for an electronically mediated interaction, 
with a polarization which is a regular function of to (see also below) . 

Now, it is easy to show analytically |t5| that in the low T limit for fi > e v n 

Im£ R {k,t) = -7rafe(e + c) , e x < e < e 2 , (5) 

where ak = So~<? ^k—q(®)> c = l l ~ e vH, and e\ = fi — e v u, £2 = £c + £vH — and e c < 8t is a characteristic energy. 
Similarly p5[ , in the high temperature limit T > (jx— c v h)/£ Jl6[ we obtain 

ImT, R (k, e) = -2ita k T . (6) 

The result is the same for /j, < c v h, save for a multiplication by 2 of the r.h.s. of eq. (5). We note that the T and 
e dependence of the result are independent of k - thus leading necessarily to a linear in T resistivity and a 
linear in 1 /e optical conductivity, even with inclusion of vertex corrections in the calculation. 

The overall behavior of ImV(q,uj) = Imx {q,u) |e(g, w)| 2 , where e(q,u>) is a properly defined dielectric function, 
follows from any screened interaction between the carriers. Hence the argument for the linear in energy and tempera- 
ture behavior of the scattering rate r _1 (T, e) is equally generic. It relies essentially on a large coefficient for the linear 
in energy term of Imxo, and the presence of vHs near the Fermi surface. What is more, in our numerical solution we 
observe that the energy e V H of the singularities is pushed by the interactions close to fi. The shape of the imaginary 
part of the self-energy imS(fc,e) of the interacting system being responsible for the modification of the density of 
states N(e), through the relation N(e) = —Tr ImG(k, e)/ir. Namely, JmS(fc, e) has a peak below /z and a dip above 
it, which account for the transfer of the spectral weight. This feedback effect reinforces the role of the singularities, 
which are more effective in producing a linear scattering rate the closer they are to the Fermi surface. This seems 
to be a plausible explanation for the common characteristic of a good many cuprates whose vHs are located between 
10-30 meV below the Fermi surface (Dp (see also Q). 

It is interesting that the electron doped Nd2- a; Ce 2 ;Cu04+5 whose vHs are at approximately ^(-350 meV [[IS), has a 
usual Fermi liquid r _1 (T) = const. T 2 [^). This lends support to the picture described above. 

Markiewicz has given a review of related work in the frame of the so-called van-Hove scenario ]2(J . 

The celebrated Millis-Monien-Pines susceptibility 0j 

XMMp{q,u) = — , (7) 

has been used to fit the low energy part of the susceptibility of the cuprates in both NMR rate and inelastic neutron 
scattering (INS) experiments. The short range antifcrromagnetic order, a remnant of the parent antiferromagnetic 
materials, with correlation length £, is responsible for the peak of the susceptibility for q near Q. Typically £ is of the 
order of the lattice constant (£ decreases as the doping increases), while u>sf ~ 10 — AQmeV. 

The origin of the small magnitude of ujsf has remained elusive thus far. E.g. Chubukov, Sachdev and Sokol 
have interpreted it as a damped spin wave mode p]]] . Spin waves are clearly well-defined excitations in underdoped 
cuprates. However to date there is no experimental proof that they do so for small energies in the normal phase of 
the optimally doped and overdoped regimes. 

An alternative explanation is a fermionic origin for LUgp. First, let us consider the total susceptibility put forward 
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Xt(q,u) = - — — . (8) 

1 - Jq{XAF(q,u) +XF{q,u)) 

Xaf(?) is due to the localized Cu spins, which interact via the effective spin exchange J q , and Xf is a purely fermionic 
susceptibility. This Xt(q,u) encompasses in an appealing way the idea of the entangled carrier-spin dynamics in the 
cuprates. Now, we use Xt above with 
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XAF{q,u) = -— -pr, 1 - 7 r^ 7TT ' XF(q,u>) = Xo(l + iuj/uj + 0(lu 2 )) , q^Q . (9) 

If |/(w)| <C w/w and 1 > ^Q$ 2 Xi we recover essentially xmmp(?,w) - which is itself an approximate form of the true 
susceptibility - with lusf — ► w cx w . 

The origin of small w (g) = Vefe F <f ||| is the proximity of the vHs to the Fermi surface - c.f. the discussion 
in the context of the scattering rate above. From the numerical solution of our system, we easily obtain values of 
lj comparable to the experimentally relevant ones, when the chemical potential is near e V H- Moreover, when e v u 
approaches fi, oj is quickly suppressed towards zero. As already mentioned, the small difference [i — e v u is found 
in several cuprates - see jlTj . It has also been suggested that this characteristic may be true irrespective of the 
doping, as long as the latter is appropriate for superconductivity. It is then clear that the carriers residing in the 
vicinity of the vHs, with a large density of states, give rise to a small uj as discussed above. Hence it is very interesting 
to know how universal this band-structure characteristic of the cuprates is, as it may explain naturally the magnitude 
of u> . This idea can be tested experimentally by e.g. determining by INS the value of u) cx ui for Nd2-zCe :r Cu04+,5, 
which is expected to be enhanced as a result of the big difference between e V H and fi. 

Moreover, we thereby propose that the effective non-interacting susceptibility of the carriers is given by the following 
ansatz 

Xo(q,u) - > Xo // (<?, w ) = Xo(q,u) + axAF(q,u) . (10) 

Xo and xaf are as above and < a < 1 is a weighting factor (which should in general depend on doping, band- 
structure etc.). The intuitive idea is that the carriers should feel the ordered antiferromagnetic background of the 
CuC>2 planes, even in the absence of phase separation j27j,^8|. In this way, the fermionic susceptibility acquires an 
antiferromagnetic enhancement, which may then influence pairing, through an electronically mediated interaction - 
see below. Note that with this effective \ e Jf our many-body approximation remains conserving, since the relation 
£ = SF/5G is valid, as xaf does not depend on G. 

Further, Xc/^ can explain the temperature dependence of the Hall resistivity of the cuprates p9| , as it can produce 
an effective carrier potential V peaked at Q. 

To locate T c we solve the gap equation 

A(fc, e n ) = V ^ k ~ k '> e " " e ^ G(k '> e «)l 2A ( fc '< e «) ■ ( n ) 



The pairing potential V p = V x (q,u m ) + V h (q,uj m ) ,V x (q,uj m ) = 1 -( Ao +u^xi(q^ rn ) i+A aXo {q^ m ) ' Vfl = 
(A +u ) 2 Xo ( q ^ m ) a 
l-(A„+I7o)xo(g,wm) B- 

The highest T c 's correspond to a d x 2_ y 2 gap, and are obtained by including an antiferromagnetically enhanced 
susceptibility in the calculation, following our ansatz eq. (10). We consider two similar forms for the AF susceptibility, 
namely (A) : XApil^m) = X a X^Tr 1 6>(w c - |w m |), where r, = £~ 2 + (q - Qi) 2 , and (B) : Xaf iq,^m) = 
x o D ELiWm ~ (2/7r)w m arctan(w m /£>) - T t D + (2/ir)TiD arctan^-)}/^ - (r.D) 2 ], with Q t = (±tt, ±tt) and 
D being a cut-off frequency, above which ImxAF(q,^) = 0. Form (B) has appeared in pq ]. Here the characteristic 
spin wave frequency los cx £~ z , and the z = 2 scaling regime has been assumed. 

With t — 2507716^, A a — and ^ = 1 (increasing the latter leads to a reduction of T c ) we have the following results. 
For form (A) we obtain T c ~ 125°K, for aX Q = 0.2eF _1 , ui c = At, t' = —0.1K, t" = 0M, U Q = 1.21 eV and n = 0.88. 
For form (B) we obtain T c ~ 105°^, for aX a = Q.beV- 1 , D = 16*, *' = -0.11*, t" = 0.25*, U D = lAleV and n = 0.91. 

We note that the basic effect of considering moderate values A a < is to allow for greater values of U before the 
AF instability condition U to tXo = 1, with U to t — A a + U , sets-in (recall that we take x%^ peaked at Qi = (±tt, ±7r)). 
Hence a finite phonon interaction allows for a bigger on-site/paramagnon coupling. 

There is an optimum value of a for the maximum attainable T c . This is due to the form of the pairing potential 
V p above : the AF instability condition allows UtotXo < 1 only. Now, for a given value of the latter product, the 
highest V p - which in principle yields the highest T c as well - will correspond to the highest possible Utot- This in turn 



corresponds to a smaller Xo, and hence to a small but finite optimum a. In fact, Zheng et al. |26| noticed that the 
highest T c 's are obtained for a combination of both optimum total carrier concentration and of a reduced relaxation 
rate (1/Tx)c u °c T Yvai u ^^2 q A 2 Imx A F(q,u) / 'w , in properly normalized units (here A q is the hyperfine coupling). 
These conclusions are in accordance with our picture, which yields both special values of the filling factor n, as well 
as special small values of the product cixaf as a prerequisite for the highest attainable T c 's. 

Monthoux and Pines have obtained similarly high transition temperatures with their approach. 

We note that the vHs at q Q suppress a s-wave gap if the pairing potential V p >0 (V p defined above is negative 
for sufficiently large and negative A a and/or U Q , and positive otherwise). The s-wave gap of Nd2- a; Ce, r Cu04+,5 can 
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be due either to V p <0 for some relevant parts of the phase space in this material, or to the fact that the vHs are 
350meV below the Fermi surface, and hence ineffective here, or possibly to both facts. Yet another possibility which 
may coexist with the above characteristics is that V p is strongly peaked close to zero momentum owing to the band 
structure of this material. All these factors can result in the s-wave T c being higher than the d-wave T c , and hence in 
the dominance of the former over the latter. 

In brief, let us discuss the possible connection to the physics of the underdoped cuprates, which are believed to be 
in a phase separated regime, with AF domains of spins separated by stripes of holes |^,^8|. One can envisage that 
with doping increasing towards the optimal regime, the stripes melt into an effective Fermi liquid, and the physics 
described here is recovered. 

To summarize, we present a single plane Fermi liquid model which for the normal state can explain the salient 
transport properties, the low energy dependence of XmmPi an d their relation to the existence of van-Hove singular- 
ities close to the Fermi surface. We also obtain enhanced transition temperatures T c . Special attention is paid to 
Nd2- a; Ce a ;Cu04+5, the properties of which, despite appearances, we believe to be fully consistent with those of the 
majority of cuprates. 

The author has enjoyed discussions with Yia-Chung Chang, Gordon Baym, Joseph Betouras, Girsh Blumberg, Nejat 
Bulut, Antonio Castro Neto, Lance Cooper, Tony Leggett, Dmitrii Maslov, Philippe Monthoux, Jorg Schmalian, Raivo 
Stern, Qimiao Si and Branko Stojkovic. This work was supported by the Research Board of the University of Illinois, 
the Office of Naval Research under N00014-90-J-1267 and NSF under DMR-89-20538. 
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